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Abstract 

Using three coupled harmonic oscillators, we present an amplitude-amplification method for 
Q-i factorization of an integer. We generalize the method in [arXiv:1007.4338j by employing non- 

orthogonal measurements on the harmonic oscillator. This method can increase the probability 
of obtaining the factors by repeatedly using the nonlinear interactions between the oscillators and 
non-orthogonal measurements. However, this approach requires an exponential amount of resources 
for implementation. Thus, this method cannot provide a speed-up over classical algorithms unless 
its limitations are resolved. 

PACS numbers: 03.67.Ac 
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I. INTRODUCTION 

Quantum computing based on qubits has attracted considerable attention (see, e.g., []]- 
Q]). There are several candidates to realize quantum computers, such as using nuclear spins 
in molecules, photons, trapped ions, superconducting circuit and quantum dots (see, e.g., 
0,0]). However, it is still a great challenge to build a large-scale quantum computer. 

Quantum computers can significantly outperform classical computers in doing some spe- 
cific tasks [3 5, 9]. For example, two important quantum algorithms are Shor's jlo| and 



Grover's [11] . Shor's algorithm 10j can factorize a large integer in polynomial time, of- 
fering an exponential speed-up over classical computation. Grover's algorithm [lj| gives a 
quadratic speed-up in searching database. This search algorithm has been found to be very 
useful in other related problems {4], 0,0. To date, the study of quantum algorithms is a 
very active area of research (see, e.g., [9j). 



Using three coupled harmonic oscillators, we have recently proposed 
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an alternative 



approach for factoring integers. We consider these three harmonic oscillators to be coupled 



together via nonlinear interactions [12] . To factorize an integer N, this approach involves 
only three steps: initialization, time evolution, and conditional measurement. In this ap- 
proach, the states of the first two harmonic oscillators are prepared in a number-state basis, 
while the state of the third oscillator is prepared in a coherent state. The states of the first 
two harmonic oscillators encode the trial factors of the number N. The nonlinear interac- 
tions between the oscillators produce coherent states that simultaneously rotate in phase 
space with different effective frequencies, which are proportional to the product of two trial 
factors [3]. In this way, all possible products of any two trial factors can be simultane- 
ously computed, and then they are "written" to the rotation frequencies of the coherent 
states in a single step. The resulting state of the first two oscillators is the factors' state 



121 ] by performing a conditional measurement of a coherent state rotating with an effective 
frequency which is proportional to N. However, the probability of obtaining this coherent 
state becomes low when N is large. In this paper, we can circumvent this limitation by 
using an iterative method for increasing the chance of finding the states of the factors. This 
amplitude-amplification method involves a number of iterations, where each iteration is very 
similar to the factoring approach we recently proposed [l^ . 

Now we briefly describe this amplitude-amplification method for factorization using three 



coupled harmonic oscillators. Let us now consider the first step of our approach. Initially, 
the first two harmonic oscillators are in a number-state basis and the third oscillator is in a 
coherent state. Let the three coupled harmonic oscillators evolve for a period of time. The 
detection is then conditioned on a coherent state with a rotation frequency being propor- 
tional to N. The probability of finding this coherent state can be adjusted by choosing both 
an appropriate period of time evolution and magnitude of the coherent state. Here we find 
that this probability is not small. Indeed, the probability of finding the factors' state can 
be increased by a factor which is the reciprocal of the probability of obtaining this coherent 
state. But the word "probability" in these sentences is different to the total probability of 
obtaining the factors. 

The resulting states of the first two oscillators, after the first step, are used as new 
input states in the second step of our approach. Also, the state of the third oscillator is 
now prepared as a coherent state with the same, or higher, magnitude. By repeating the 
same procedure described in the first step, we can obtain the states of the factors with a 
much higher probability. We then iterate these procedures L ~ (log 2 N) times, until the 
probability of finding the factors' state is close to one. 

As an example of how this method works, we show how to factorize the integer 1, 030, 189 = 
1009x1021. Here the probabilities of obtaining coherent states, with rotation frequencies 
proportional to N, are larger than 0.1 in each iteration. The probability of finding the fac- 
tors can reach nearly one after 12 iterations. By comparing with the examples N = 101, 617 
and 10,961, we can show that the required number of iterations for factoring logarithmically 
scales with N. 

However, this approach requires an exponential amount of resources for its implementa- 
tion. First, it requires an exponential amount of energy to encode a number N onto the 
state of a harmonic oscillator. The input energy scales with the number N. Therefore, 
the required energy becomes enormous if the number N is large. Second, it is necessary to 
prepare an exponential size O(N) of ensemble of oscillators for this iterative approach. This 
is because a large number of oscillators are abandoned after a conditional measurement. A 
number at least of order of O(iV) oscillators are thus required to be prepared for the input 
size N. Therefore, this approach does not provide any speed-up for factorization compared 
to classical algorithms. 

This paper is organized as follows: In section II, we introduce a system of coupled har- 
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monic oscillators. In section III, we study the quantum dynamics of the coupled harmonic 
oscillators starting with a product state of number states and a coherent state. In sec- 
tion IV, we propose an amplitude-amplification method to factorize an integer using three 
coupled harmonic oscillators. We discuss the convergence and performance of this factoring 
algorithm. For example, we show how to factor the number 1, 030, 189 using this approach. 
In section V, we discuss the problems and limitations of this approach. Finally, we make a 
summary in section VI. 



II. SYSTEM 

We consider a system of (A + B) coupled harmonic oscillators. The Hamiltonian of the 
j-th harmonic oscillator is written as 

H 3 = —2- + 3 3 3 m 

where Uj is the frequency of the harmonic oscillator, rrij is the mass of the particle, and 
j = 1, . . . , A+B. The operators Xj and Pj are the position and momentum operators, which 
satisfy the commutation relation, [Xj,Pj] = ih. The annihilation and creation operators of 
the j-th harmonic oscillator are defined as 



a , p^i y. + ^L I (2) 



4 = M(x>-£r)- (3) 



2h \ rrijUJ 



The operators aj and aj obey the commutator [a^aj-] = 1. The Hamiltonian of the three 
harmonic oscillators can be expressed in terms of the annihilation and creation operators: 

A+B 

H = h 2J Uja^jCij, (4) 
i=i 

Here we have ignored the constant term. 



We consider the harmonic oscillators coupled to each other via nonlinear interactions 
Such nonlinear interactions can be described by the Hamiltonian Hj as 



12|. 



A+B 

Hi = h ^ fk{ni, n A )a\a k (5) 

k=A+l 



where fk(n>x, . . . ,ua) are linear or nonlinear-operator functions (excluding divisions) of the 
number operators rij = a^aj, for j — 1, . . . , A, and k — A + 1, . . . , A + B. 
The total Hamiltonian H = H + Hj can be written as 

A+B A+B 

H = u i a ) a i + ^ ^ fc ( ni ' • • • ' n A)a[a k . (6) 

j=l k=A+l 

The Hamiltonians H and Hj commute with each other, i.e, [H , Hj] = 0. The total Hamil- 
tonian H in Eq. (jH]) is exactly solvable. The eigenstate \E mi mA B ) of the Hamiltonian H 
is a product state of number states \rrij)j of the harmonic oscillators, i.e., 

A+B 

\E mi ,...,m A+B ) = IJ l m i)i> ( 7 ) 
i=i 

corresponding to an eigenvalue £ mii ... iBA+s 

= w i m i + J] ^( mi > • • • ' m A) m k- (8) 
i=i fc=A+i 

III. QUANTUM DYNAMICS IN PHASE SPACE 

We study the time evolution of the fc-th harmonic oscillator in phase space starting 
with a state FT. \ m j)j ® \ a )ki where FJ . |mj)j is the product state of the number states of 
the harmonic oscillators, and \a) k is the coherent state of the k-ih harmonic oscillator, for 
j = 1, . . . , A, and A; is a number from A+l to A+B. By applying the time-evolution operator 
U{t) = exp(-iHt) [H is the Hamiltonian in Eq. §6§] to the initial state FT. \ m j)j\ a )k, it 
becomes 



U(t) Y[ \m j ) j \a)^ 



exp 



^Ujm^jt Y[ ">.i) 



3 3 



x exp { - i[u k + / fe (mi, . . . , m A )}a\a k t} \a) k , 



exp 



j j 
where a m i,...,m A (^) is a complex function, i.e. 



- i {^2 UJ i m 3^j t Yl\ m j)j\ a mi! ..., mA (t))k, (9) 



am u ..., mA (t) = exp[-iQ mi! _ imA t]a, (10) 




FIG. 1. (Color online) Schematic diagram of the time evolution in phase space of the coherent 
state of the harmonic oscillator k, corresponding to the product of number states FJ ■ \mj)j, where 
j = 1, . . . , A. The coherent state with complex amplitude a is here depicted as a light blue circle in 
phase space (Xk,Pk), where X k and P k represent the position and the momentum of the harmonic 
oscillator k. After a time t, the coherent state rotates about the origin with an angle t£l mi ... mA . 



and Q 



mi,...,mA 



is an effective rotation frequency of the coherent state in phase space 



n 



m\,...,mA 



w fc + /& (mi, . . .,m A ). 



(11) 



In Eq. ©, we have used the relation 



cxp 



(-i-da\a k )\a) k = |aexp(-ii?)) fc , 



(12) 



where d is a phase factor. 

Note that the product state of the number states Ylj \ m j)j ls an invariant; namely it 
does not change with time. Nonlinear interactions, described by the Hamiltonian Hi in 
Eq. (jSJ), cause the coherent state of oscillator k to rotate about the origin with a frequency 



n 



mi,...,ra^ 



in phase space. From Eq. lull) , the frequency fl 



L m\,...,mj\ 



depends on the number 



states Ylj \ m j)j °f the A oscillators. A schematic diagram of the time evolution of the 
coherent state of the harmonic oscillator k in phase space is shown in Fig. [TJ 



IV. FACTORIZATION 



Extending our proposal in Ref. 12j, we now present an amplitude-amplification method 
to factor any positive integer N. For example, let us consider three coupled harmonic 
oscillators for factorization by setting A = 2 and B = 1. Using this method, the probability 
of finding the factors' state can reach nearly one after L iterations. Simplified schematic 
diagrams of the factorization approach are shown in Figs. [2] and [3j 



Amplitude Amplification 



Measurement 



FIG. 2. (Color online) Factoring algorithm using three coupled harmonic oscillators. Initially, the 
state p^ of the harmonic oscillators 1 and 2, is prepared in a number-state basis. By repeatedly 
applying L times the iterations I\ (as shown in Fig. [3]), the resulting state of the harmonic oscillators 
1 and 2 becomes the state of the factors of N. Finally, the factors of the number N can be obtained 
by measuring the state of the oscillators 1 and 2. 



Ii iteration 



\a {l) )3 3(a (l) \ 



UAti) 



Conditional 
Measurement 



(*i))3»<a»(t,) 



FIG. 3. (Color online) The Z-th iteration, J/, in factorization. In each iteration, the third harmonic 
oscillator is prepared in a coherent state |a^))3 with a magnitude |o^|. Then, a unitary operator 
Uf(ti) is applied. By performing a conditional measurement of the coherent state \orS(ti))3, the 
reduced density matrix pi 1 can be obtained. 
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A. Algorithm 



Now we consider the total Hamiltonian Hf of the system 

3 

H f = h^uja^aj + fiffaaor (ni, n 2 )4a 3 , (13) 

where the nonlinear-operator function /factor { n i, n 2) for computing the product of any two 
trial factors is written as, 

K 

/factor ("i, n 2 ) = ^2 9k (nin 2 ) k . (14) 

k=l 

Here the parameters are the coupling strengths of the nonlinear interactions (nin 2 ) k : c4a 3 , 
for k = 1, . . . , K. This operator function /f ac tor will output an eigenvalue which is a power 
series of the product n x m, for the product state \n,m) = \n)\ \m)2- 

The total state of the first two harmonic oscillators, in the number state basis, is initially 
prepared, i.e., 



E Pn™'\n,m)(n',m'\, (15) 



P 

n,n'm,m' 



where the p™™ are the probabilities of the states of the first and second oscillators, while 
n, n', m, ml — 2, . . . , N/2. The states of the harmonic oscillators 1 and 2 can be prepared in 



arbitrary states, including pure states or mixed states 12| . The states of oscillators 1 and 2 
encode all trial factors of the number N. Each number state \n)j represents each trial factor 
n, for j = 1, 2. 



1. First iteration: I\ 

In the first iteration, the state of the third harmonic oscillator is prepared in a coherent 
state la^ 1 ^. By applying the time-evolution operator U f(t{) = exp(—iH fti) to the initial 
state p (1) (0) = p ( r ] g> |«W) 3 3 (a (1) |, it becomes 

pW(*0 = UfitJpW ® | a W) 8 .(a^lujih), (16) 
= E ^V,^)l«i 1 i l (^i))3 3 («i:L,(ti)|(n',m'| (17) 

n,m,n',m' 

where 

jtfffa) = exp{z[(n' - n ) Ul + (m> - m)^}^™' . (18) 
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Here a4m(i), for I — 1, in Eq. flT7|) is a complex function, 

= exp(-in nm t)a (l \ (19) 
and fi nm is an effective rotation frequency of the coherent state in phase space 

K 

tt nm = u 3 + 9k(nm) k . (20) 

k=l 

Note that the product of the two factors of r and s is equal to N: r x s = N. The rotation 
frequency Qn of the coherent state is 

K 

Q N = u 3 + J2g k N k . (21) 

k=l 

If the product of any two numbers is not equal to N, then the frequencies f2 nm are different 
to the frequency Ojv- Thus, the state of the harmonic oscillator 3 can act as a "marker" for 



the states of factors and non-factors 



12|. 



Now we define a measurement operator Mi which can be written as [Tj 

Mi pP)(t,) = JiEfiffliti) = Y,Ei pWfa) El (22) 

i 

Ei = |a#(f0>8 3<a#fe)|, (23) 

where a$(ti) = exp(— zOjvi/) and / is the number of iterations. A conditional measure- 
ment M\ is performed on the third oscillator at the time t\. The probability of obtaining 
this coherent state |a^(ti)) 3 becomes 

Pi(E 1 ) = TT[M lP ^(t 1 )}, (24) 

E„™ \M | 2 (oft) 

n,m 

where the coefficient e™ = 3(0$ (ti)\oinm(ti))s, for I — 1, and 

e^ = exp{-|««| 2 [l-exp(^^-^ nm ^)]}, for / ^ 1, (26) 

is the overlap between the two coherent states |a^(t;)) 3 and \anm(ti)) 3 , respectively. Note 
that the value of the probability Pr(£' 1 ) can be adjusted by appropriately choosing the 
evolution time and the magnitude | ct^ 1 ^ | . In practice, this probability Px(E{) cannot be 
adjusted to be extremely small. 



Recall that a state is conditioned when this state is conditional on the measurement of a 



certain state 



17| . After the measurement of \a N (ti))s, the density matrix of the conditioned 



state can be written as 



3: 



M) = (27) 

= £■ E K^i)^^K^)(<^ , i®i«!i ) (ti))3 3(«!; ) (ti)i, (28) 

n,m,n ,m 

where C\ is a normalization constant 

w i / v fnm I nml ' \ ) 

n,m 

Note that the trace of this density matrix Tr[pW] is equal to one. After the first iteration, 
the probability of finding the factors of the number N is increased by a factor Cj~ , which 
is the inverse of the probability Pr(Ei) as seen from Eqs. ( |25i) and ( |28l) . The probability 
amplification [see also Eqs. fl53l) and (15^1) ] is thus inversely proportional to the probability 
of obtaining the coherent state ~Px(E\). 



2. Second iteration: I 2 

(2) 

After the first iteration, we now obtain the reduced density matrix p r of the first two 
harmonic oscillators as 

P? ] = ^ E pS 1 '(ti)effi l eS&|n. ( 30 ) 

n,m,n,m' 

We consider the state in Eq. f l3"Ul) of the oscillators 1 and 2 as an input state for the second 
iteration. The coherent state of the third harmonic oscillator is prepared in a coherent state 
\a^)s, with a magnitude \a^\. The nonlinear interactions between the three harmonic 
oscillators are then turned on for a time £2- The state evolves as 

P (2) (t 2 ) = ^ £ K^d^^ (31) 

n,m,n',m' 

where £2 = h + £2- 

Next, a conditional measurement .M 2 is applied to the system at the time t 2 . The 

(2) 

probability of obtaining the coherent state \a N (t2))3 becomes 

Pr(£ 2 ) = Tr[AV 2) (t 2 )], (32) 

= J_\^„nm I (1) |2 I (2) |2 /oo\ 
£i / y rnw rnml rnml " 

^ n,m 
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The conditioned state can be written as 
(2) = M 2P V(t 2 ) 

Pc Tr[M 2P W(t 2 )Y W 
= ^ E J^'&)e&^^ (35) 

n,m,n',m' 

where C2 is the normalization constant, 

°2 — / ^Pnml £ reml I Sim I ■ 



n,m 

i-l 



The probability of finding the factors is enhanced by a factor C 2 after the second iteration. 

The coefficients \&hm\ is less than one for any product n x m ^ N. From Eqs. (1291) and 
f[3T)j) . it can be seen that C2 < C\. Therefore, the probability of finding the factors is now 
higher after one additional iteration. 



3. L-th iteration: II 

Similarly, we now iterate the procedure (L — l) times. After (L — l) iterations, the reduced 
density matrix of the oscillators 1 and 2 can be written as 



pl L) = ^— E Pfiffc-O ffeSSL $* m >\n,m){n',m'\, (37) 



where 



n.m.n'.rn' 1=1 



L-1 



Cl-1 - E^™» II l e nm| 2 ' ( 38 ) 



n,m 1=1 
L-1 



Z=l 



The state of the third harmonic oscillator is now prepared in a coherent state |a^ L ^3, with 
a magnitude |a^ L ^|. Let the three coupled harmonic oscillators evolve for a time t^. This 



gives 

L-1 



P {L) (t L ) = -^- E 11^ ^m'l»^>l«S2(*i)>3 s<«^(*i)|«m'|. (40) 



n,m,n' ,m 1=1 



By performing a conditional measurement |a^(t.r,))3, the state becomes 

= (41) 
Pc Tr[.M L pW(t L )]' l4ij 

= 7^ E ^'(^JIsls^^ (42) 

n,m,n',m' 1=1 
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After the L-th step, the probability of finding the factors is increased by a factor C L X . From 
Eq. (1381) . the probability of obtaining the coherent state |o!^(£l))3 can be written as 

1 L 

^{E L ) = —Y,PnZj[\^L\ 2 , (43) 

n,m (=1 

= (44) 

The entire iterative procedure is now completed. The convergence and performance of 
this method will be discussed in the following subsections. 



B. Convergence 

We now study the convergence of this iterative method. We first consider the magnitude 
of the coherent state \a^\ for each iteration as 

\a^\^\a^\^...^\a^\^...^\a^\. (45) 

Thus, we have 

1 > |e (1) I 2 ^ |e (2) I 2 ^ . . . ^ |e (/) I 2 ^ . . . > |e (L) | 2 . (46) 

For any product of n and m being not equal to N, the coefficient |enm| 2 is less than one 



and decreasing for higher I, and the evolution time ti is non-zero and appropriately chosen. 
When the number of iterations L tends to infinity, the product of the coefficients |e 
tends to zero 



(0 12 
nm 



'i 

L 



limni^r^O. (47) 



L— >oo 

1=1 

The coefficients \er}\ 2 are equal to one for any product of two factors r and s being equal to 
N, i.e., when r x s = N. 

Now we consider the probability of finding a pair of factors, r and s, after the l-th 
iteration, which is 

Pr<(M» = ^p. (48) 

Since the coefficient |enm| 2 is less than one, the normalization constant C\ in Eq. ( )38l) is 
decreasing, i.e., G\ < Cj_i. Therefore, the probability of finding a pair of factors increases 
after an additional iteration. This shows that this iterative method is convergent. 
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From Eqs. (!38|) and ( HTj) . it is very easy to show that 

L 



— *on L— S>00 z — ' - LJ - 

n,m i=l 

= ( 5 °) 

In the limit of large number of iterations L, we can obtain the state p/ of the factors 

lim -> P/ , (51) 

where 

P/ = ^T E ?f/M<TVI, (52) 

/ r,s,r',s' 

and = XlrsPrs' with r,s,r',s' are factors of N (rxs = r'xs' = N). This shows that 
the state of the factors can be achieved by employing this iterative method, if a sufficiently 
large number of iterations is used. 

C. Performance 

We can now estimate the number L of the iterations required to achieve a probability of 
order of one for factoring N. We investigate the amplification ratio A; of the two probabilities 
of finding the factors Pr/_i(|r, s)) and Pr/(|r, s)J after the (I — l)-th and the l-th iterations. 
From Eq. ( HH]) . we have 

Note that this ratio A/ is just the reciprocal Pi -1 (Ei) in Eq. ( )44l) of the probability of 
obtaining the coherent state. Practically, this probability Pr(Ei) cannot be too small. For 
example, let the amplification ratio A; be roughly equal to A ~ O(10 _1 ) for each iteration, 
and let the probability of the factors' state Pr (|r, s)) before the iterations be of order 
0(N~ Z ), where z is a positive number. After L iterations, the probability of finding the 
states of the factors can be increased by a factor \ L . Here we require 

0{X L ) ~ 0{N Z ). (55) 
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Therefore, we obtain that the number L of necessary iterations is of order of 0(log 2 N). In 
the limit of large L, the probabilities Pt(E{) in Eq. (jHJ) tend to one because Cl approaches 
the sum of probabilities of the factors' states in Eq. (I5"0"j) . The probability of finding the 
factors will be slowly increased after the number L ~ 0(log 2 N) of iterations is reached. 



D. Example: Initial pure states 

In this section, we study how to factorize an integer iV with an initial pure state using 
this factoring algorithm. We consider the initial state of the first two harmonic oscillators 
as the superposition of number states, i.e., 



*(«)> = 7T Y l»>i TT Y hM. (56) 




where 

D 1 = {\VN]-2) 1 ^, (57) 



D 2 = ( fiV/3] - \VNTl] + 1) 1/2 , (58) 

are two normalization constants. Here we consider trial factors from 3, . . ., N/3. The 
probability of finding the product of two factors is of order of 0(N~ 3 ^ 2 ). 

For example, now we show how to factor the integer N — 1, 030, 189 = 1009x1021. For 
simplicity, we now take K = 1, which is the lowest order of nonlinearity. The Hamiltonian 
can be written as 

3 

Hi = h^^cUjOjaj + hga\aia\a2a\as, (59) 
i=i 

where g is the nonlinear strength. The stronger nonlinear strengths and high-order nonlin- 
earity can significantly shorten the required time evolution of the system 12j. But the role 
of nonlinearity is not directly relevant to the number L of the required iterations for the 
amplitude amplification. 

We take t* as the evolution time for the Z-th iteration, 

t*i = —r h (60) 
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TABLE I. This table shows the fidelities Fi, the probabilities Pr(Ei) of obtaining the coherent 
states |a^(^))3 in Eq. (f23j) . and the evolution time t* for the I iterations. Here E\ is the density 
matrix of the coherent state \a^J (t*))s in Eq. (|23p and is measured in units of g^ 1 . 



Iterations 

7 

I 


b ldelities 
1*1 


Probabilities for coherent states 


Evolution times 

H 


1 


2.010x10 


0.143 


1.704 


2 


1. 403x10 


0.143 


1.342 


3 


9.782x10 


0.143 


5.000 


4 


£? ooi * , ~i r\ — fi 

6.821x10 


0.143 


4.610 


5 


4.739x10 


0.144 


r\ too 

0.732 


6 


o orn, , 1 n — 4 

3.259x10 


0.145 


3.108 


1-7 

( 


2.172x10 


0.150 


1.635 


8 


1.445x10 


0.150 


4.559 


9 


1.045xl0 _1 


0.138 


4.222 


10 


5.092X10" 1 


0.205 


6.046 


11 


8.506X10- 1 


0.599 


2.434 


12 


9.919X10" 1 


0.858 


1.175 


13 


9.985X10- 1 


0.994 


5.089 


14 


9.997X10" 1 


0.999 


5.833 


15 


l.OOOxlO" 1 


1.000 


0.708 



where r> is a uniformly distributed random number on the interval [0, 1]. We now evaluate 
the performance of this method by investi gati ng the fidelity F\ between the reduced density 



matrix pr +1 ^ and the factor's states pj as 



18 



IS | 



Tr 



{pfp^pf) 1 ' 2 



(61) 



Table [I] shows the relevant fidelity, the probability of obtaining the coherent state 
\o${t*))% with a magnitude \a®\ = 2, and the evolution time t* for each iteration. Before 
starting the iteration, the initial fidelity Fq is very low: 2.883 x 10~ 9 . In the first few steps, 
the probabilities, Pr(i^), of obtaining the coherent states, are about 10" 1 . We emphasize 
that now the probabilities Pt(Ei) are not extremely small, even when N is large. This 
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FIG. 4. (Color online) Bar chart: the average of the fidelities between the reduced density matrix 
Pr and the factors' states pj are plotted versus the number I of iterations in (a) N = 1,030, 189, 
(b) N = 101,617 and (c) N = 10,961. In each iteration, the fidelity F\ is taken at the time ft, 
which is a uniformly distributed random number ranging from to 2n/g. Here the sample size 
is equal to 100 fidelities. The bars, which indicate the mean values, are shown in light blue. The 
error bars, which indicate the standard deviations, are shown in black. The dashed lines indicate 
that the fidelities attain 0.9. 



resolves the limitation of our previous proposal 12) . The probability of finding the factor's 
state can be increased by 10 after a single iteration. After ten iterations, the fidelity can ex- 
ceed 0.5. The fidelity can reach nearly one after two more iterations, while the probabilities 
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P(Ei) approach one. Also, we have explicitly shown that this method can factor a number 
of order of 10 6 with 12 iterations. 

Since the evolution time t* in Eq. (160 p is a uniformly distributed random number from to 
27r/g, it is necessary to study the effect of this randomly chosen time t* to the performance 
of factoring. Now we study the averages of the fidelities Fi and probabilities Pt(Ei) of 
obtaining the coherent states. We also examine their standard deviations in each iteration. 

In Fig. S(a), we plot the average of the fidelities Fi versus the number I of iterations. Here 
we take the sample size to be equal to 100 fidelities in each iteration. As shown in Fig. HJ the 
average of the fidelities is about 0.5 at the 10-th iteration. The fidelity is greater than 0.9 
at the 12-th iteration. After 12 iterations, the mean values of the fidelities reach nearly one. 
The error bars, which indicate the standard deviations, are shown in the same figure. We 
can see that, after the 12-th iteration, the standard deviations are relatively small compared 
to the mean values in each iteration. 

We then plot the average of the fidelities Fi versus the number I of iterations for iV = 
10,1617 and 10,961 in Fig. @Ja) and (b), respectively. The fidelities exceed 0.9 after the 
10-th iteration for N = 101,617 and the 8-th iteration for N = 10,961. The numerical 
results in Fig. [H show that the required number of iterations increases logarithmically with 
TV for which the fidelity is greater than 0.9. 

In Fig. [3(a), the average of the probabilities Pr(i^) of obtaining the coherent states is 
plotted versus the number I of iterations. In the first ten iterations, the probabilities Pt(Ei) 
are about 0.1. Then, it increases and saturates around one after 13 iterations. Moreover, 
these standard deviations are much smaller than the mean values of Pt(Ei). This means 
that the statistical effect of the time t* t is small on the performance of quantum factorization. 

In Fig. [5(b) and (c), the average of the probabilities Pr(Ei) of obtaining the coherent states 
are plotted versus the number I of iterations for N = 101, 617 and N = 10, 961, respectively. 
The probabilities are greater than 0.9 after the 11-th iteration for N = 101,617 and the 
9-th iteration for iV = 10, 961. Therefore, the results show that the number of iterations 
(Pr(^)>0.9) also logarithmically scales with N. 
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FIG. 5. (Color online) Bar chart: the average of the probabilities of obtaining the coherent states 
\or^f(t*))s are plotted versus the number I of iterations in (a) A?" = 1,030,189, (b) A" = 101,617 
and (c) A?" = 10, 961. In each iteration, the probability is taken at the time tl which is a uniformly 
distributed random number from to 2ir/g. Here the sample size is taken to be 100 probabilities. 
The bars, which indicate the mean values, are in red. The error bars, which indicate the standard 
deviations, are shown in black. The dashed lines indicate that the probabilities attain 0.9. 

V. LIMITATIONS AND PROBLEMS 



We have introduced an amplitude-amplification method for factoring integers by using 
three coupled harmonic oscillators. However, this method requires an exponential amount 
of resources for implementation. This means that the required resource is of order of 0(N), 
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where N is the input size. This means that this approach cannot provide any speed-up com- 
pared to classical algorithms. We discuss these two limitations in the following subsections. 

A. Exponential energy resource 

First, this approach requires an exponential amount of energy resource to encode a num- 
ber N onto the state of a single harmonic oscillator. For example, to factorize a number 
N, two harmonic oscillators are used for encoding 0(iV 3//2 ) possible states of trial factors. 
The energy E ~ HuO(N 3 ^ 2 ) is thus required, where u ui 1:2 are the frequencies of the 
harmonic oscillators 1 and 2. The required energy scales with the size of the input number 
N. Therefore, the energy becomes enormous when the encoded number is large. This en- 
coding method becomes impractical when a large input number is used. This problem may 
be resolved by using another encoding method. For example, the number could be encoded 
either onto qubit or qudit states. Thus, the required energy for encoding numbers could be 
reduced. 

B. Exponential size of the ensemble 

Second, this method requires an exponential size of the ensemble of harmonic oscillators 
for an input N. In each iteration, it is necessary to abandon a large number of harmonic 
oscillators after each conditional measurement. The number of abandoned oscillators is 
proportional to the failure probability of obtaining conditional measurement of the coherent 
states in each step. Therefore, a number at least of order of 0(N) harmonic oscillators are 
needed in order to complete the entire procedure. This requires an exponential resource 
for preparing the ensemble of oscillators. This approach may be improved by employing 
efficient methods for preparing the ensembles. 

VI. CONCLUSIONS 

We have presented an amplitude-amplification method by repeatedly using the nonlinear 
interactions between the harmonic oscillators and non-orthogonal measurements. We have 
shown that this approach can be used for factoring integer, and the factors of an integer N 
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can be obtained, with a high probability, by using a number of iterations L ~ 0(log 2 A). 
We have numerically studied an example for factoring iV = 1,030,189, respectively. We 
have shown how to factorize an integer of order of O(10 6 ) within 12 iterations. In each 
iteration, the probability of obtaining this coherent state, with the rotation frequency being 
proportional to N, is not less than 0.1. By comparing examples with iV = 101,617 and 
N = 10, 961, we have shown that the required number of iterations increases logarithmically 
with N. 

However, using coupled harmonic oscillators, this method requires the use of an expo- 
nential amount of resources, i.e., exponential energy and ensemble size. Thus, this approach 
becomes impractical for large input sizes. We hope that the resolutions could be proposed 
to overcome the problems of this approach. 

Also, we stress that the nonlinear interactions between the coupled oscillators and condi- 
tional measurements are essential in this approach. By appropriately controlling nonlinear 
interactions between the coupled harmonic oscillators, the functions with integer inputs can 
be evaluated in a single operation. To implement this approach, it is necessary to engineer 
"many-body" interactions of the system of harmonic oscillators. For example, to perform 
quantum factorization, it is required to generate "three-body" interactions between the har- 
monic oscillators. We have briefly discussed the possible implementations in Ref. 12]. One 
of the promising candidates is neutral atoms or polar molecules trapped in optical lattices 



20j-|22j. The "three-body" interactions can be tuned by external fields 20|, |21 |. 
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